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Recent progress in denning jets 

Gregory Soyez a * 

a Brookhaven National Laboratory, Building 510, Upton NY 11973, USA 



From dedicated QCD studies to new physics background estimation, jets will be everywhere at the LHC. 
In these proceedings, we discuss two important recent series of improvements. In the first one, we introduce 
new algorithms and new implementations of previously existing algorithms, in order to cure limitations of their 
predecessors and to satisfy fundamental requirements. In the second part, we show that it is of prime importance 
to carefully choose the jet definition — algorithm and parameters — to optimise kinematic reconstructions at the 
LHC. Noticeably, we show that while at scales around 100 GeV, R ~ 0.5 is an appropriate choice, clustering at 
the TeV scale requires R ~ 1 for optimal efficiency. We finally show that our results are valid in the presence of 
pileup, provided that a subtraction procedure is applied. 
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1. Introduction 

In collider physics, as soon as the final state in- 
volves hadronic particles, jets become fundamen- 
tal objects present in many studies. Even if the 
idea of a jet as a bunch of collimated partons or 
hadrons is what one always keeps in mind, the 
concept of a parton itself is ambiguous and as a 
consequence different jet definitions exist, possi- 
bly leading to different sets of jets for the same 
hadronic event. 

In practice, many jet definitions, which we shall 
review in Section [21 have been used for recent 
analysis. It turns out however that, among them, 
some fail to satisfy the fundamental requirements 
agreed upon in 1990 (Section [272]) . 

The first series of results presented in these 
proceedings addresses those failures and give so- 
lutions that have been proposed recently. This 
includes the FastJet implementation of the k t al- 
gorithm, Section [2731 as well as the the SISCone 
and anti-fc* algorithms, Sections 12.41 and 12.51 Fi- 
nally, in Section l2~6l we will briefly discuss a fil- 
tering technique using jet substructure that has 
been recently proposed to reduce sensitivity to 
the underlying event. 

Since one has the choice between different jet 
definitions to perform jet analysis, a legitimate 
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question is which of them is best suited for a 
given analysis one wants to perform. In Section 
[3] of these proceedings, we answer that question 
in the case of simple kinematic reconstructions at 
the LHC. To quantify the efficiency of a jet defini- 
tion, we introduce a figure of merit that is directly 
related to an effective luminosity ratio. We will 
see that not being careful enough in the choice of 
the jet algorithms and its parameters (basically 
the radius R) can lead to important consequences 
for potential discoveries at the LHC. 



2. Meeting fundamental requirements 

2.1. Defining jets in the 20 th century 

To begin with, let us briefly review the algo- 
rithms that have been widely used over the past 
two decades for jet reconstruction in pp collisions. 
Generally speaking, they fall in two categories 
that we discuss hereafter. 

Successive recombinations. The first fam- 
ily of jet clustering algorithms works by defin- 
ing a distance between any pair of objects and 
a beam distance for every object. One identifies 
the smallest distance; if it is a beam distance, the 
object is called a jet and removed from the event, 
otherwise, the two objects are recombined in a 
single one. The procedure is repeated until no ob- 
ject are left in the event. Two well-known exam- 
ples of recombination algorithms are the kt |H2j 



1 



2 



Gregory Soyez 



and Cambridge/ Aachen (C/A) 3 algorithms, us- 
ing the distance 

da = min(fc t 2 ^fc^)(A4+A4), 

where p = 1 (p = 0) corresponds to the fcj (C/A) 
case. 

Cone. The cone algorithms aim at defining 
jets as dominant directions of energy flow. To 
achieve that goal, one defines the concept of stable 
cone as a circle of fixed radius in the (y, <fi) plane 
such that the sum of the 4-momenta of the parti- 
cles inside it points in the direction of its centre. 
Most of the cone algorithms used so far are seeded 
in the sense that the search for stable cones starts 
from a given set of seeds from which one iterates 
the cone contents until it is stable. Since stable 
cones might overlap, one cannot define them as 
jets directly. Two different techniques are in use 
to overcome this overlapping problem. 

The first option, that we shall refer to as cone 
algorithms with split-merge, first identify the set 
of all stable cones, then run a split-merge proce- 
dure on them. The latter repeatedly identifies the 
two hardest (originally, in E t ) overlapping cones; 
if their overlap passes a threshold fraction (the 
overlap parameter /) of the softer of the selected 
cones, they are merged. Otherwise, they are split 
by associating every particle to the cone to whose 
centre it is closer. 

For the seeded versions of the cone algorithms 
with split-merge, one usually starts with all the 
particles in the event (usually with a pt threshold) 
as seeds, as is the case for the CDF JetClu [4] and 
ATLAS Cone algorithms. A step forward is to 
add as new seeds the midpoints between all pairs 
of stable cones found after this first pass. This is 
the case for the majority of the recently-used cone 
algorithms [5], noticeably the CDF MidPoint, DO 
run II Cone and PxCone algorithms. 

The second solution to the problem of overlap- 
ping stable cones, iterative cone algorithms with 
progressive removal (IC-PR), starts by iterating 
the stable cone search from the hardest seed in 
the event. The resulting stable cone is called a jet 
and its contents are removed from the event. One 
then proceeds by iterating from the hardest re- 



maining seed, until all the particles arc clustered. 
The characteristic feature of this type of cone al- 
gorithm is that it produces hard jets that are cir- 
cular and soft-resilient. In other words, the addi- 
tion of soft particles does not modify the shape of 
the hard jets, which is sometimes seen as an ad- 
vantage for calibrating the jets. The CMS Itera- 
tive Cone algorithm [6] falls in this sub-category. 

2.2. Fundamental requirements 

Back in 1990, a list of fundamental require- 
ments that every jet definition has to fulfil was 
agreed upon [7J. It is known as the SNOWMASS 
accords and consists of the 5 following criteria: 

1. simple to implement in an experimental 
analysis; 

2. simple to implement in the theoretical cal- 
culations; 

3. defined at any order of perturbation theory; 

4. yields finite cross section at any order of 
perturbation theory; 

5. yields a cross section that is relatively in- 
sensitive to hadronisation. 

As consequences, because of constraints [T] and 
[5l we want the implementation of the algorithm 
to be fast enough and as insensitive as possible to 
the underlying event (UE) so that it can be used 
in experimental analysis. Also, the requirement 
that the cross section remains finite at any order 
of perturbation theory implies that the algorithm 
has to be infrared and collinear (IRC) safe. In- 
deed, if it was not the case, cancellation between 
real emissions and virtual corrections would not 
happen properly, leading to divergences in per- 
turbative cross-sections. 

2.3. Speeding up the k t algorithm 

Beside its rather large sensitivity to the UE, 
one argument that was sometimes used against 
the k t algorithm was its relatively slow running 
time. The KtJet 8 and KtClus pQ implementa- 
tions, of complexity O (A" 3 ) where N is the num- 
ber of particles in the event, have a clustering 
time around 1 second for N — 1000, compared to 
about 0.2 seconds for the MidPoint cone with a 1 
GeV seed threshold. 

Remarkably, the complexity can be reduced 
to O (Nlag(N)) using computational-geometry 
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technique^] [S]. As seen on figure [T] (k t (Fast Jet) 
curve), this leads to a considerable improvement. 

2.4. SISCone and the IR unsafety of the 
MidPoint algorithm 

It had already been noticed that the JetClu 
algorithm had some IR unsafety problems e.g. 
when two hard particles were distant by more 
that R and less that 2R, an additional infinitely 
soft gluon added between the particles could 
change the clustering from 2 jets to 1 jet, yielding 
unreliable cross-sections at NLO in the inclusive 
jet cross-section. The MidPoint algorithm was 
then introduced to cure that problem. 

Unfortunately, this is not the end of the story. 
Even if situations with 2 hard particles in a com- 
mon vicinity (plus another one, hadronic or elec- 
troweak to balance pt) are IR safe provided one 
uses the MidPoint algorithm instead of JetClu, 
the problem has just been shifted to situations 
with 3 hard particles in a common vicinity. This 
is illustrated in figure [21 where we have clustered 
a 3-hard-particle event twice — with, fig. E^b), 
and without, fig. [Ha) an additional infinitely soft 
gluon — and notice that one finds different sta- 

2 For N < 5000, there exists a O (N 2 ) implementation 
that turns out to be faster than the TVlog(iV) [9]. 



Figure 2. Stable cones found by the MidPoint 
algorithm for (a) a 3-particle event, (b) the same 
event with an additional infinitely soft gluon. 



ble cones in the two cases. This means that the 
MidPoint algorithm is also IR unsafe, though one 
order further in the perturbative expansion in the 
strong coupling than the JetClu algorithm. 

To solve that problem to all orders in the per- 
turbative expansion, we first notice that the IR 
unsafety comes from the fact that, in the event 
without the soft gluon, the stable cone enclos- 
ing particles 2 and 3 has been missed. Since the 
mathematically well-defined set of stable cones is 
IRC safe, i.e. it changes neither when splitting 
a particle collinearily nor when adding infinitely 
soft particles (up to harmless stable cones made 
only of those soft particles), finding a method 
that provably identifies all stable cones, guaran- 
tees IRC safetjU. 

3 We also need to be careful about the IRC safety of the 
split-merge procedure. Actually, a couple of technical is- 
sues, such as the choice of the variable used to order the 
overlapping cones, have to be dealt with. We will not 
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We then observe that every circular enclosure 
of given radius R in the (y, (f) plane can be trans- 
lated in any direction until it touches one point, 
then rotated around that point until it touches a 
second one, without changing its contents. There- 
fore, enumerating all pairs of points, and for ev- 
ery pair of points considering the two circles of 
radius R they define and the four possible in- 
clusion/exclusion states of the edge particles, we 
enumerate all possible enclosures. For each of 
them we can then test if it is stable or not, which 
solves our problem. The complexity of this al- 
gorithm is O (TV 3 ) (a factor N 2 coming from the 
enumeration of the pairs of parent points and an 
additional factor of N to test the stability of ev- 
ery enclosure). It is actually possible, using ex- 
tra geometric observations, to improve that com- 
plexity to O (N 2 log(iV)) . This has been imple- 
mented jlOj in a new algorithm named SISCone 
(Seedless Infrared Safe Cone). Fig. Q] shows that 
SISCone runs faster than the IR-unsafc MidPoint 
algorithm (of order iV 3 ), even when a seed thresh- 
old of 1 GeV is applied. 

2.5. Anti-fc f and the collinear unsafety of 
the IC-PR algorithm 

In a similar way as for the case of the Mid- 
Point algorithm, we can show that the iterative 
cone algorithm with progressive removal (IC-PR) 
also suffers from divergences in the perturbative 
series, at the same order as MidPoint, this time 
due to collinear unsafety. In order to see that, 
first consider the event of fig. 02a). Iterating from 
the hardest seed gives one jet containing all parti- 
cles. If one splits the hardest of these particles in 
two collinear ones, fig. (3^b), iteration starts with 
the leftmost particle and 2 jets are found. This 
means that the IC-PR is collinear unsafe at the 
level of 3 particles (+1 to balance pt). 

To address this issue, we will go back to the 
recombination- type algorithms. We have already 
mentioned that setting p = or 1 in eq. |T]) re- 
produces the C/A and k t algorithms respectively. 
We now introduce the anti-k t algorithm, corre- 
sponding to p = — 1 in eq. JT]) [TT] . 

At first sight, it is not obvious what this new, 
IRC-safe, algorithm has to do with the IC-PR. 
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Figure 3. Jets found by the iterative cone for 
(a) a 3-particle event, (b) the same event with a 
collinear splitting. 



However, since hard particles will be associated 
a small anti-fct distance, they will grow in circles, 
clustering softer particles in their vicinity up to a 
distance R. This thus leads to soft-resilient hard 
jets, i.e. the boundary of the hard jets is not af- 
fected by soft radiation, the precise characteristic 
feature of the IC-PR. 

Finally, the anti-fct algorithm is also amenable 
to a fast implementation (see fig. [Q) using the 
same techniques as for the kt algorithm (see Sec- 
tion O). 



address them here (see 1101 for details). 



2.6. Filtering 

One of the promising potential improvement of 
jet clustering is to make use of the jet substruc- 
ture. For example, we can use the following fil- 
tering technique to reduce the contamination due 
to the underlying event: 

1. Cluster the event using, e.g., the C/A algo- 
rithm with a radius R, 

2. For each jet, recluster it using a smaller ra- 
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dius i? su b and keep only the n su b hardest 
jets as part of the initial jet, throwing the 
other subjets. 
The aim of filtering is to remove contamination 
due to soft background like the underlying event 
while keeping as much as possible of the pertur- 
bative radiation. This has already proven [12j to 
be efficient in Higgs searches in the bb channel. 

In what follows we shall use -R S ub = R/2 and 
rt su b = 2, though a more extensive study of the 
effects of these parameters would be interesting. 

3. Quantifying kinematic reconstruction 
efficiency 

Now that we dispose of 5 IRC-safe jet algo- 
rithms — the k t , C/A, anti-fc t , SISCone and 
C/ A+filtering algorithms, all available from Fast- 
Jet [13] — we may ask, given an analysis involv- 
ing jets we want to perform, which jet definition, 
i.e. the jet algorithm and its parameters, is best 
suited. In this Section we address that question 
for the case of kinematic reconstructions at the 
LHC. 

We first introduce a series of benchmark pro- 
cesses we will investigate, then a figure of merit 
that allows one to quantify the performance of a 
jet definition and finally present our results, both 
with and without including pileup. For a more 
extensive discussion, see |14|15j . 

3.1. Benchmark processes 

We will study the following 3 processes: 

• Z' — > qq as a source of quark jets. We re- 
construct the Z' from the 2 hardest jets in 
the event (imposing a maximal rapidity dif- 
ference | Ay | < 2 between them). 

• H — > gg as a source of gluon jets. We re- 
construct the H from the 2 hardest jets in 
the event (imposing again |Ay| < 2). 

• a ti pair with fully hadronic decay into 6 
jets (4 jets from the 2 W bosons and 2 jets 
from the b and b quarks). The two Ws and 
then the 2 tops are reconstructed from the 
6 hardest jets in the evenlQ. 

4 Practically, 6-jets are tagged assuming the b mesons are 
stable; then the two W are reconstructed by pairing the 
4 remaining jets so as to minimise (Mi 1 i 2 — My/) 2 + 



For the case of quark and gluon jets, we can study 
the scale dependence by varying the mass of the 
Z 1 / H boson (in practice, between 100 GeV and 
4 TeV). The last case, fully hadronic ti decay, 
allows us to study the relevance of our results 
in more complex environments, where one risks 
tensions between resolving the jets and capturing 
the perturbative radiation. In each of these situ- 
ations we have generated a sample of events using 
Pythia 6.4 tune DWT. 

3.2. Figure of merit 

In order to quantify which jet definition is per- 
forming better than another, we need a figure of 
merit. Since jets are expected to represent an 
original parton, one might be tempted to quan- 
tify the efficiency by comparing the jets to the 
initial partons. But because partons are a ill- 
defined concept (especially at NLO), this is not 
robust enough. 

Another option would be to fit a given distri- 
bution, e.g. a Gaussian, to the peak in the mass 
spectrum. But since the shape of the peak can 
be asymmetric, this would not produce a reliable 
result cither. 

For these reasons, we will quantify the effi- 
ciency of a jet definition based on the observation 
that, if two peaks have similar number of events 
then the narrower is the better one. We therefore 
introduce the quality measure QJ =Z as the width 
of the smallest possible window that contains a 
fraction f — z of the events. The smallest val- 
ues of QJ =Z should then correspond to the most 
efficient jet definitions. 

This intuitively does what we want in the sense 
that a "better" jet definition should contain a 
given fraction of the events in a smaller window 
and therefore have a smaller QJ =Z - 

This quality measure can be related to a varia- 
tion of luminosity needed to maintain a constant 
significance for a signal relative to background. 
For a jet definition JD, the latter is defined as 

E(JD) ee A^^v^bkgd- Assuming a constant 



(Mj„j 4 — Myy) 2 ; finally, the two top jets are reconstructed 
by matching the fo-jets with the W so as to minimise the 
mass difference between the 2 top candidates. 



6 



Gregory Soyez 




Figure 4. result for the quality measure for quark jets at 100 GeV (left) and gluon jets at 2 TeV (right). 
The results are displayed as a function of R and each curve corresponds to a given algorithm. 



background, we thus have 

1/2 
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two jet definitions JDi and JD 2 . This relation 
comes from the fact that the number of signal 
events is fixed by the fraction / = z in the quality 
measure, while the background is directly propor- 
tional to the width of the window i.e. to QJ= Z - 
A "better" definition, i.e. a smaller QJ =Z , thus 
corresponds to a larger discriminating power. 
We can then define an effective luminosity ratio 



p £ (JD 2 /JD! 



£ (needed with JD 2 
£(needed with JDij 

E(JDi) f Q?U(JD2) 
S(JD 2 )J Q7 =z (JDi) 



(2) 



This means that a jet definition JDi with a qual- 
ity measure twice as large as JD 2 will need twice 
the integrated luminosity in order to achieve the 
same discriminating power as JD 2 . 

3.3. Results without pileup 

For the processes presented in Section [3T| we 
have clustered our event samples using the 5 IRC 
safe jet algorithms and varying the parameteJl 
R between 0.1 and 1.5. We can then compute 



5 In the case of SISCone, the overlap threshold / has been 
fixed to 0.75, the preferred choice at the time being. 



the quality measure in each of those cases (see 
figure SI), which allows us to (i) compare different 
algorithms, and (ii), for a given algorithm, find 
the optimal value i?best of the parameter R. 

Figure [5] summarises our results in a more com- 
pact and physical way: we have plotted the effec- 
tive luminosity ratios corresponding to a series of 
selected processes: quark and gluon jets at 100 
GeV and 2 TeV, and top reconstruction in the 
tt sample. For every process, pc has been nor- 
malised to the best jet definition for that process. 

We observe a few important features: first of 
all, in general, SISCone and C/A with filtering 
tend to perform slightly better than the kt, C/A 
and anti-fct algorithms, especially at higher scales 
(with the exception of the top reconstruction 
where all algorithms have similar performances). 
For example, for the 2 TeV gluon case, choosing 
the kt algorithm instead of SISCone or C/A with 
filtering which are the preferred choices, trans- 
lates into a cost of nearly 50% in the effective 
luminosity ratio. 

Furthermore, the efficiency of a jet definition 
strongly depends on R and not choosing the pre- 
ferred value i?best can be very costly. On top of 
that, i?bcst varies significantly from one process 
to another. It increases from i?bcst — 0.5 at small 
scales up to i?bost — 1 at TeV scales. i?best is also 
bigger for gluon jets than for quark jets. Practi- 
cally, using SISCone with R = 0.5, the preferred 
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Figure 5. Effective luminosity ratios as a function of the parameter R. Each line corresponds to a different 
process and each column to a different algorithm. For each process, pc is normalised to the best possible 
definition for that process. 
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choice for quark jets at 100 GeV, to cluster gluon 
jets at 2 TeV costs a factor of 2 in pn compared to 
the best definition for that process. Conversely, 
using the preferred value for gluon jets at 2 TeV 
- i.e. SISCone with ft = 1.1 or C/A with filter- 
ing and ft = 1.3 — to cluster quark jets at 100 
GeV also leads to pc increasing by at least 50%. 

3.4. Results with pileup 

Pileup corresponds to the fact that multiple pp 
interactions can happen at the same time. For the 
case of the LHC at the designed luminosity, one 
typically has an average of about 25 collisions per 
bunch crossing. This produces a large number of 
additional soft particles that form a reasonably 
uniform background over the detector. 

For our concerns, pileup has two consequences. 
Because it adds a background to every jets, the 
pt of the jets will be overestimated and the po- 
sition of the mass peak in our kinematic recon- 
structions will be shifted towards larger masses. 
Also, because the amount of pileup varies from 
one event to another, the reconstructed peak will 
be smeared, an effect that directly affects our 
quality measure. 

We therefore want a (simple and generic 
enough) method to subtract the contamination 
from the pileup background on an event-by-event 
basis. In practice, we will follow the subtraction 
method suggested in |16j : for each jet j in an 
event, we first compute its (4-vector) area [17] 
Aj. The subtracted jet is then obtained using 

p£ 8 „b=p?-M?, (3) 

where p is the average density of pileup per unit 
area. Since the density of background is reason- 
ably uniform, we estimate p for each event using 
16 p = median{ptj/j4 tJ }, where all the jets (up 
to a maximal rapidity) are included in the com- 
putation of the median. 

In practice, we have considered the same 
benchmark processes as in the study without 
pileup. Pileup is added to each event under the 
form of a random (Poisonian) number of mini- 
mum bias events, also generated with Pythia 6.4 
tune DWT. We apply the same reconstruction 
procedure as above and compute the quality mea- 
sure QJ =Z with pileup effects subtracted or not. 



We illustrate our results by showing on figure 
[6] the effective luminosity ratio pc as obtained for 
a few representative processes. We compare the 
results with pileup (subtracted or not) with the 
corresponding quality measures obtained without 
pileup (solid curves). 

The first observation is that if the pileup is 
not subtracted (dashed curves), it causes a large 
degradation of the quality measure. The pre- 
ferred ft is also shifted to smaller values, which 
one can explain by the fact that the contamina- 
tion due to pileup background is smaller at small 
ft. The second element of information comes from 
the case where we apply our subtraction method 
(dotted curves). Obviously, even if the subtrac- 
tion is not perfect — the quality is still larger than 
before pileup addition — one sees a significant 
improvement compared to the situation without 
subtraction. This means that the subtraction 
method gives narrower peaks i.e. the smear- 
ing due to pileup fluctuations between events is 
strongly reduced by the subtraction. On top of 
that, after the subtraction has been performed, 
the preferred value for ft obtained in the situation 
without pileup is no longer strongly disfavoured 
to the profit of a smaller ft. This is important 
since it implies that our conclusions from Section 
13.31 are still valid in the presence of pileup pro- 
vided one uses subtraction. 

4. Conclusions 

The first point we have addressed concerns the 
jet algorithms themselves. As summarised in ta- 
ble [T] some of the commonly used algorithms fail 
to satisfy the fundamental requirements of in- 
frared and collinear safety We have introduced 
SISCone and the anti-fct algorithm to cure in- 
frared unsafeties of the MidPoint algorithm and 
collinear unsafeties of the IC-PR. We have also 
mentioned a new implementation of the k t algo- 
rithm leading to a sizeable improvement in speed. 

We have then addressed the problem of quan- 
tifying the performances of the jet definitions in 
the case of kinematic reconstructions at the LHC. 
From that study, one learns important messages 
showing that optimisation of the jet definition can 
significantly improve the potential for an early 
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Figure 6. Same as for figure now including the effects of pileup. The solid (red) curve corresponds to 
the situation without pileup, the dashed (green) one to the case of unsubtracted high-luminosity pileup 
and the dotted (blue) curve to the situation where pileup is subtracted. 
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Algorithm 


Type 


IRC 


inclusive kt [1121 


SR p= i 


OK 


Cambridge/ Aachen [3] 


SR p= o 


OK 


anti-fc* UJJ 


SR p= _i 


OK 


SISConc [10, 


SC-SM 


OK 


CDF JetClu 0] 


IC r -SM 


IR2+1 


CDF MidPoint [5J 


IC„ ip -SM 


IR-3+1 


DO Run II cone [5] 


IC mp -SM 


IR-3+1 


ATLAS Cone 


IC-SM 


IR2+1 


CMS Iterative Cone [6] 


IC-PR 


Coll 3+ i 



Table 1 

Overview of some jet algorithms used recently in ex- 
perimental or theoretical work. SR p=:r = sequential 
recombination (with p — 0, ±1, see |T])); SC = seed- 
less cone (finds all cones); IC = iterative cone (with 
midpoints mp, ratcheting r), using either split-merge 
(SM) or progressive removal (PR) in order to deal 
with overlapping stable cones. Regarding IRC status, 
given n hard particles in a common neighbourhood, 
IRn+i indicates that the addition of 1 extra soft par- 
ticle can modify the number of final hard jets, while 
Coll n +i means that the collinear splitting of one of 
the particles can modify the number of final hard jets. 



discovery at the LHC. We have shown that the 
preferred jet definition varies significantly from 
one process to another. SISCone and C/A with 
filtering tend to perform slightly better than the 
other algorithms and the preferred value for R 
goes from R ~ 0.5 for 100 GeV jets to R ~ 1 
for TeV jets with a slightly larger value for gluon 
jets compared to quark jets. This indicates that 
a single choice for R is not sufficient to cover the 
whole kinematic range at the LHC. 

Finally, if one includes the effects of pileup, 
we have seen that, provided we use an adequate 
subtraction procedure, the conclusions obtained 
in the analysis without pileup hold once pileup 
is added, noticeably, there is no need to take a 
smaller value of R. 

In the more complex case of top reconstruc- 
tion where jet clustering involves tension between 
catching enough radiation to reconstruct the jet 
energy, and resolving the 3 jets coming from the 
top decay products, techniques like the subjet 
analyses might be of crucial importance and def- 
initely deserve further studies. 
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